The Bg and Dg decay constants in 3 flavor lattice QCD 
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Capitalizing on recent advances in lattice QCD, we present a calculation of the leptonic decay 
constants /s^ and fo^ that includes effects of one strange sea quark and two light sea quarks. 
The discretization errors of improved staggered fermion actions are small enough to simulate with 
3 dynamical flavors on lattices with spacings around 0.1 fm using present computer resources. By 
shedding the quenched approximation and the associated lattice scale ambiguity, lattice QCD greatly 
increases its predictive power. NRQCD is used to simulate heavy quarks with masses between 1.5mc 
and nib. We arrive at the following resuhs: fs, = 260 ± 7 ± 26 ± 8 ± 5 MeV and fo, = 290 ± 
20 ± 29 ± 29 ± 6 MeV. The first quoted error is the statistical uncertainty, and the rest estimate the 
sizes of higher order terms neglected in this calculation. All of these uncertainties are systematically 
improvable by including another order in the weak coupling expansion, the nonrelativistic expansion, 
or the Symanzik improvement program. 



PACS numbers: 12.38.Gc, 13.20.Fc, 13.20.He 



We present the first complete calculation of the Bg and 
Dg decay constants with 3 flavors of sea quarks with small 
masses. The Ds decay constant, /d, , has been measured 
P, 0] , and CLEO-c promises to reduce the experimental 
errors significantly. The comparison of experimental and 
lattice results will be a vital test of lattice QCD. The Bs 
and Bd decay constants, fs, and /s, are necessary in 
order to constrain Vtd via — B^ mixing. Since experi- 
mental measurements of Jb, remain elusive, an accurate 
lattice calculation is critical for improving phcnomeno- 
logical tests of CKM unitarity. 

The calculation presented below makes use of lattice 
Monte Carlo simulations, done by the MILC collabora- 
tion which include the proper sea quark content: one 
dynamical strange quark and two flavors of dynamical 
quarks with masses as light as ms/4. The correct num- 
ber of flavors is necessary in order for the lattice theory 
to have the same /3- function as real QCD. Only then can 
we expect, in principle, lattice results to agree with ex- 
perimental measurements. 

Inclusion of light up and down sea quarks is essen- 
tial for accurate lattice phenomenology. The innova- 
tion which allows this on present computers is an im- 
proved staggered discretization of the light quark ac- 
tion H H 1, S H S [IS For each dynami- 
cal flavor the fourth root of the fermion determinant 
is used in the Monte Carlo updating algorithm. De- 
spite some open theoretical issues concerning staggered 
fermion algorithms, these calculations are free of am- 
biguities present in quenched simulations, provided the 
u and d sea quark masses are light enough. In prac- 
tice, calculations involving unstable particles or multi- 
hadronic final states require new or refined techniques 
before they can be reliably simulated; however, there 
are several quantities for which numerical calculation is 



straightforward. Ref. [iJl presents lattice results for a va- 
riety of these "gold-plated" quantities, including /jr and 
Jk, and several splittings in the T spectrum. The change 
from Uf = Q results, which differ substantially from ex- 
periment, to n/ = 3 results, which agree with experi- 
ment up to the 3% lattice uncertainties, suggests that the 
improved staggered fermion method correctly simulates 
QCD. Heavy-light pseudoscalar leptonic decay constants, 
presented here, fit into the "gold-plated" category 12]. 

Our calculation of Jb^ and fu, uses standard lattice 
QCD methods. Correlation functions were computed us- 
ing a subset of gauge field configurations generated by 
the MILC collaboration. The configurations include the 
effects of 2 dynamical light quarks with equal bare mass, 
denoted by rrvf^, and 1 dynamical strange quark, with 
bare mass mf^. The lattices we used have spacings of 
about 1/8 fm and volumes of about (2.5 fm)^ x 8.0 fm. 
We focus on the configurations where mf'^/rrvf''^ — 1/5 
and 2/5. As we discuss below, the physical strange quark 
mass, rris, obtained from the light hadron spectrum is 
actually 4/5 of nvf^, so the light sea quark masses are 
approximately rrig/A and mgjl. 

Details of the Monte Carlo simulations which gener- 
ated the ensemble of gauge fields are given in In 
this work we use a level splitting in the T spectrum, e.g. 
the T(2S — IS) splitting, to determine the lattice spac- 
ing, instead of the length scale t\ , derived from the static 
quark potential, which was used in Q. The T spectrum 
will be presented in detail in a future publication [l3 |: 
however, we note here that as the light sea quark mass 
is increased, the T(1P — IS) splitting becomes slightly 
smaller than experiment when the lattice spacing is set 
by T(2S - IS). With statistical errors between 1 - 2%, 
the difference between experiment and lattice values for 
the IP — IS splitting is Oct, Icr, and 2a for rwf^ /rris ~ 
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TABLE I: Heavy quark parameters and meson masses for two 
values of the light sea quark mass. 



aMo 


n 


mf « ms/4 (568 config) 
m^^ (GeV) m^"f (GeV) 


^ ms/2 (468 config) 
m^^ (GeV) m^";* (GeV) 


2.8 


2 


5.6(2) 


5.3(4) 


5.22(17) 


5.4(4) 


2.1 


4 


4.38(10) 


4.2(3) 


4.28(11) 


4.3(3) 


1.6 


4 


3.52(6) 


3.5(3) 


3.56(7) 


3.5(3) 


1.2 


6 


2.84(5) 


2.78(19) 


2.93(4) 


2.83(19) 


1.0 


6 


2.53(4) 


2.41(16) 


2.60(3) 


2.45(16) 



1/4, 1/2, and 3/4, respectively. Except for estimating 
sea quark mass effects, we use tlie mf^'^/ms « 1/4 lattice 
to obtain our results. 

Tlie analysis of Q has also been updated with re- 
spect to determining the quark mass corresponding to 
the physical strange quark mass sector. Rather than us- 
ing ss mesons which are either unstable (the 0) or do not 
exist (the pseudoscalar) we use the results of the partially 
quenched chiral perturbation theory analysis of uik and 
in El- On lattices with auomf^ = 0.007,0.01, 
and 0.02, the physical strange sector is obtained with 
awom™' = 0.040. (Here we keep explicit a factor of the 
mean-link uo, which was absorbed into the definition of 
the hght bare quark mass in j^j.) 

Tabled lists the bare heavy quark masses used in this 
work, along with the corresponding NRQCD stabiliza- 
tion parameter n. For each mass, we estimate the heavy- 
light meson mass two ways: the "kinetic mass" m^" is 
extracted from finite momentum correlators using the 
meson dispersion relation, and the "perturbative mass" 
m^"^* is estimated from the zero momentum correlator 
and the one-loop heavy quark mass renormalization. We 
find that the bare heavy quark mass uMq = 2.8 produces 
the correct experimental mass for the T 0| and the Bs 
within statistical errors. 

When we construct correlation functions for heavy- 
light mesons, we use the equivalence between staggered 
and naive fermion propagators. The operators we use 
are the same ones as with Wilson-like discretizations or 
in the continuum limit. This method for computing the 
Bs mass and decay constant has been presented, along 
with tests on quenched lattices, in recent work |l6l| . 

The Bs decay constant is defined through the axial 
vector matrix element; (0|A^|i3s(pp)) = fs^Pn- Here we 
use only the temporal component. Up through 0{1/Mq) 
three lattice operators contribute to Aq, the leading-order 
current, Jq, and two sub-leading currents Ji and J2; the 
matching from the lattice to the continuum is done in 
perturbation theory ^ one-loop calculation yields 

Ao = (l + a.p„)j(")+(l+a,pi)j(^'^"^'+a,p2J^''^"'' (1) 
where jj^^'^^^^ = Jq^-* — asCioj!f \ and so too for J^^'^^^K 




0.2 0.4 0.6 

l/m„ (GeV-i) 



FIG. 1: Squares show = fus y/'mRs vs. I/ttihs on the 
nfi'^ /nis ~ 1/4 lattice. The solid line shows the fit described 
in the text, and the asterisk shows the value of extrap- 
olated to the charm sector. Experimental values for l/m_B^ 
and 1/mDs are shown as dashed vertical lines. 



These operators explicitly subtract, through one-loop, 
the power-law mixing of Jq°^ with Jp^^ and Jq^^ 0|. The 
perturbative calculation which determines poi Pi, P2, CiOi 
and C20 will be presented separately [l^. 

From fits to correlation functions, we extract the com- 
bination (^H, = fH^y/rfiH^. Let us denote by the 
contribution of J^^ to ^Hs- For this calculation the 
matrix elements are computed for mesons at rest, so 
$(1) ^ $(2)^ Table in summarizes fits to the numerical 
data, converted to physical units using 1/a = 1.59(2) 
GeV for the mf^ « TOs/4 lattice and 1/a = 1.61(2) 
GeV for the mf"^ « ms/2 lattice. (The quoted un- 
certainties come from statistical and fitting uncertain- 
ties in the T(2S — IS) splitting 0|.) By comparing 
(j)(i)/$(o) .^vith $(i.sub)^(j)(o) Qjjg observe the sizable 

power law mixing of J*-^^ with j'^' . The expression (jT]) 
absorbs the mixing back into the term proportional to 
J*^°\ so <i)(i'^"*^)/<i)(o) represents the physical contribu- 
tion of l/Mg terms to to two-loop corrections. 
The 4% contribution from the 1 /Mq operator we see for 
the Bs {aMo = 2.8) is the same size seen in quenched 
studies over a range of lattice spacings . 

Figure ^ shows the heavy quark mass dependence of 
on mH (plotted as squares). The data are fit well 

by' 

= *r (2) 

with a correlated per degree-of-freedom of 0.7/2. We 
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TABLE II: Simulation results for ^^f^ = /j^^ ^/niH^ for each 
light sea quark mass and heavy quark mass. The third col- 
umn lists the leading order term ^^"^ = (1 + asPo)^^°\ The 
fourth and fifth columns show the contributions of Jq'\ with 
respect to before and after the power law subtraction. 

The sixth column gives the result for ■ Statistical and 
fitting uncertainties are quoted in parentheses (not including 
the statistical uncertainty in 1/a). 









«l/4 




aMo 


-l'^0(GeV3''2) 




(^(l,sub) ^<j,(0) 


$H,(GeV='''2) 


2.8 


0.640(11) 


-9.0(4) % 


-3.7(4) % 


0.614(13) 


2.1 


0.598(10) 


-11.7(4)% 


-5.0(4) % 


0.565(11) 


1.6 


0.557(8) 


-14.7(4) % 


-6.4(4) % 


0.516(9) 


1.2 


0.519(7) 


-18.3(4)% 


-7.8(4) % 


0.470(8) 


1.0 


0.499(6) 


-20.7(4) % 


-8.6(4) % 


0.445(7) 






sea / 

m£ /ma ' 


=.1/2 




aMo 


$^0(GeV='/2) 




^^(i,sub)/^(0) $^^(GeV3/2) 


2.8 


0.640(15) 


-8.9(6) % 


-3.6(6) % 


0.615(14) 


2.1 


0.599(11) 


-11.4(6)% 


-4.7(6) % 


0.567(12) 


1.6 


0.563(7) 


-14.2(5)% 


-5.9(5) % 


0.523(9) 


1.2 


0.528(6) 


-17.7(5)% 


-7.1(5)% 


0.481(6) 


1.0 


0.510(6) 


-20.0(5) % 


-7.9(5) % 


0.459(6) 



find $5^^* = 0.85(4) GeV3/2, Ci = -1.82(20) GeV, 
and C2 = 1.59(35) GeV^. It is important to note that 
most of the mass dependence of comes through 

the action and not through the I/Mq currents. A fit 
to the leading order = (1 + a,po)*^°^ (Table ^ 

yields similar fit parameters: ^^tat.LO ^ o.83(4) GeV^/^, 
C^o ^ -1.53(19) GeV, and = 1.32(32) GeV^ 

Since the action is responsible for most of the mu^ de- 
pendence and is accurate through 0(AQ(~.j-,/mQ), per- 
forming the operator matching through 0{\/M^) should 
yield only small corrections to the values for Ci, 
and C2 listed above. 

We use the fit above to interpolate '^h, slightly to the 
physical value of l/ms^. The result is fs, — 260(7) MeV, 
where the quoted statistical error combines statistical er- 
rors from 1/a and a^^^^B^- Uncertainties in m^" are 
small compared to the other uncertainties in the fit. The 
systematic uncertainties due to the neglect of higher or- 
der terms are estimated by assuming coefficients of 0(1). 
For example, two-loop terms omitted in ^ are estimated 
to be 10% effects, taking the coupling constant defined 
through the plaquette Us = a^I {2/a) = 0.32 This 
is the largest systematic uncertainty. Leading discretiza- 
tion errors are 0{asa'^ ^^QJ-|), where Aqcd is the typical 
scale of nonperturbative dynamics. Taking 1/a — 1.6 
GeV and Aqcd = 400 MeV impHes 2% cutoff effects. 
The operator matching neglects terms 0{A^Qj^/'mD or 
approximately 1%. We note that the coefficient of the 
spin-dependent term in the action is the tree level one. 



0.3 
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^ aMo = 2.8 : 




i ^ : 

aMo=1.0 1 


1 1 1 1 
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FIG. 2: Dependence of <l>£r^ on the light sea quark mass. The 
upper data points correspond to the Bs {aMo = 2.8) and 
the lower points come from our lightest heavy quark mass, 
aMo ~ LO. The different symbols indicate which quan- 
tity was used to set the lattice spacing: octagons use the 
T(2S — IS) splitting and diamonds use the T(1P — IS) split- 
ting. Error bars represent combined statistical uncertainties 
of 1/a and o^/^^w,. 



giving rise to an 0(q:s Aqcd /toq) error in the action, 
which is about 3% for bottom quarks. This does not 
necessarily translate into a comparably large error in de- 
cay constants where the spin dependence is not obvious; 
however, to be conservative, we quote an overall 3% un- 
certainty due to relativistic corrections. 

Extrapolating the fit described above to the physical 
value for l/rriD^ requires care since higher order terms 
in l/niH^ become increasingly important. The fit above 
extrapolates to Jd^ — 290(10) MeV. We use a Bayesian 
analysis to estimate possible effects due to higher order 
terms. Allowing terms like C„/to^ with n > 3 in the 
fit (I2Jl, with Gaussian priors for C„/GeV" = zt S and 
(5 = 1 — 4, we find higher order terms could lead to a 20 
MeV error in fjj^ . This is taken as the statistical and fit- 
ting uncertainty. The 0{a^) perturbative error and the 
©(asa^Aggj-,) discretization error are again estimated 
to be 10% and 2%, respectively. For charmed mesons 
Aqcd/"^c = 1/4, so the ^(Aqqq/toq) corrections to the 
matching are estimated to be 6%. As mentioned above, 
the O (as Aqcd /'Tig) error in the spin dependent term in 
the action probably does not lead to a proportional error 
in the decay constant. Nevertheless, we quote 10% as a 
conservative estimate of possible relativistic corrections. 

Figure El shows vs. light sea quark mass for two 
heavy quark masses. We show the results where the lat- 
tice spacing is set from the T(2S — IS) splitting and the 
same data adjusted as if the spacing were set from the 
T(1P — IS) splitting. No sea quark mass dependence is 
observed for uMq = 2.8, nor is a significant dependence 
observed for qMq = 1.0 when the spacing is set with 
T(2S — IS). Only because there is a 2% ambiguity in 
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setting the scale with mf^'^ « ms/2 does there appear to 
be some sea quark niass effect in the aMg = 1.0 data with 
the T(1P — IS) spacing. No such scale setting ambiguity 
exists for the data with mf'^ « nig/ 4,^ which is where our 
main results are obtained. 

Experimental measurements of the leptonic decay con- 
stants for Bg and Dg are challenging, so these are rare 
cases where lattice QCD can lead experiment. The 
most recent and precise experimental results for Jd^ are 
280±17±25±34 MeV |l| and 285±19±40 MeV which 
agree well with the calculation presented here. Shrinking 
the experimental uncertainties on fjj^ is a major goal of 
the CLEO-c program. 

Turning to comparison with the existing literature on 
lattice QCD, we quote recent world averages and com- 
ment on another new result. The Uf — Z results pre- 
sented here have central values larger than lattice calcu- 
lations with Uf = or 2. For example, recent averages 
im of quenched results are /^f"" = 200(20) MeV and 
= 230(14) MeV, significantly lower than our re- 
suhs. The 2 flavor world averages, /^^"^ = 230(30) MeV 

and — 250(30) MeV, are higher than the quenched 

and agree within the quoted uncertainties. 

Recently the JLQCD collaboration reported a lattice 
result using two dynamical flavors of improved Wilson 
fermions with mass between 0.7tos and 2.9tos ,22]. They 
quote /s, = 215(9) (+2) (13) (^o) MeV, with the first er- 
ror statistical, second from chiral extrapolation of the 
sea quark mass, third from finite lattice spacing com- 
bined with truncation of NRQCD and perturbative ex- 
pansions, and the fourth from ambiguity in setting the 
strange quark mass. This calculation uses much larger 
quark masses than ours, uses unstable hadron masses to 
set the lattice spacing and to^, and does not include a 
dynamical strange quark. Further work will be required 
to determine which, if any, of these accounts for the dif- 
ferences with our results. For instance, it would be inter- 
esting to calculate from several T splittings on the 
configurations of '2^1 and compare /^^ using those scales 
to the value quoted in their paper, which uses the scale 
set by vfip. 

Recent sum rule calculations agree within errors for 
the Bg decay constant: e.g. = 236(30) MeV ;23;, 24]. 
On the other hand, they do not calculate an increase in 
the decay constant as the heavy quark mass decreases: 
e.g. = 235(24) MeV H. 

To summarize, we have completed a calculation of the 
Bg and Dg decay constants using 3 flavor lattice QCD. 
The more realistic sea quark content allows a unique lat- 
tice spacing to be determined using one of several quan- 
tities 4.3] , enabling more reliable prediction of quantities 
not yet measured experimentally. Our final results are 

fs, 260 ± 7 ± 26 ± 8 ± 5 MeV 

fo, = 290 ± 20 ± 29 ± 29 ± 6 MeV . (3) 



The uncertainties quoted are respectively due to statistics 
and fitting, perturbation theory, relativistic corrections, 
and discretization effects. The result for the Dg decay 
constant agrees with experimental determinations, and 
the result for the Bg decay constant is a prediction for 
future experiments. Improvement of the lattice results 
requires a two-loop perturbative matching calculation or 
the use of fully nonperturbative methods. On the other 
hand, much of the perturbative uncertainty cancels in the 
ratio Jb^ / fsd ■ Work is underway to study /b^ using the 
methods discussed in this paper [isf . 

Simulations were performed at NERSC. We thank the 
MILC collaboration for their gauge field configurations. 
This work was supported in part by the DOE, NSF, 
PPARC, and the EU. 
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